Nonequilibrium quasiparticle relaxation in superconducting qubits 
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In a superconducting qubit the lifetime of qubit state is restricted by nonequilibrium quasipar- 
ticle tunneling. We calculate the rate of these tunnelings using the nonequilibrium effects they 
induce on the condensate chemical potential of leads and islands. We show that by decreasing 
temperature below a crossover, the quasiparticle relaxation rate changes from exponential to much 
slower suppression and saturates to a finite value at zero temperature. This prediction is consistent 
with recent experiments. Our model also indicates a striking modification to qubit transitions: the 
matrix element of an energy transition in qubit strongly depends on coupling between tunneling 
quasiparticles and the environment. This features important fabrication hints to improve quantum 
state efficiency. 
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Superconducting integrated devices are among the 
leading candidate for quantum information processing 
with meticulous state control p] [2] ■ Their primary chal- 
lenge for quantum computation is their fragile quantum 
states that make them vulnerable to various types of 
decoherence mechanisms. The nature of some of these 
mechanisms are yet to be understood [3H2] • Experiments 
on superconducting quantum circuits have shown the en- 
ergy relaxation time reaches a saturation at about 1-10 
/is and may not be further enhanced by lowering temper- 
ature [BHH]- At low temperature although equilibrium 
quasiparticles are depleted yet a density of nonequilib- 
rium quasiparticles (NQP) may tunnel in and out from 
superconductors. The lack of coherent time enhancement 
is a consequence of tunneling of these excess quasiparti- 
cles. They make the qubit energy levels and correspond- 
ing states vary stochastically in time, thus degrades com- 
putational performance [91413J . 

Experiments have shown the density of excess quasi- 
particles depends on fabrication and observed in the 
range of 0.1-100//im 3 [HtHo] . Martinis et.al. in Ref. Q2] 
and 118) formulated this phenomenon by Fermi's golden 
rule in a special case of large capacity junction and nu- 
merically determined the tunneling rate in the frame- 
work of nonequilibrium thermodynamics. In Refs. [1 11 — 
113] the qubit part of the tunneling rate was generalized to 
cover junctions of arbitrary macroscopic phase; however, 
the tunnelings are restricted only to lowest-lying energy 
states. 

An equilibrium analysis |12j shows that by decreasing 
temperature the rate of quasiparticles tunneling expo- 
nentially suppresses such that a few milliKelvins cool- 
ing at low temperature results in an effective blockade. 
In experiment such a blockade has not been seen, in- 
stead the relaxation rate was seen to slowly vary at low 
temperature |15H17| . This indicates that the quasiparti- 
cle tunnelling process can have a nonequilibrium nature. 
In Ref. [19] we extended the theory treat the coupling 
between the quasiparticle channel and the environment 
nonperturbatively, hence permitting strong coupling be- 



tween these two subsystems. If, as is the case for qubits, 
the environment largely consists of the LC-mode of the 
junction, this corresponds to multi-quantum transitions 
between its eigenstates. 

Upon entry of a quasiparticle the parity of the super- 
conducting state is altered from even (in the lack of un- 
paired particle) to odd [2"UH22"] . This change of parity 
bring the condensate out of equilibrium, where the en- 
ergy distribution is no longer the Fermi-Dirac function. 
In order to determine the energy distribution one ap- 
proach is to take advantage of the quadratic dependence 
of recombination rate on quasiparticle density. A numer- 
ical analysis shows the tunneling rate calculated from 
such an energy distribution captures the experimental 
features of temperature dependence of the rate, however 
the derivation requires numerical efforts [17] , 

Here we consider the charge imbalance effect [2"3Tf2"5] 
that a superposition of quasi-electrons and holes, when 
tunnels, may induce on a superconducting region. Us- 
ing these nonequilibrium effects, we calculate the tun- 
neling rate. Our model shows that by cooling down the 
superconductor the tunneling rate from an exponential 
suppression turns into a much slower variation and ap- 
proaches a finite value at zero temperature. This be- 
havior has been observed for the rate in recent experi- 
ments |15H17j . Our nonequilibrium analysis determines 
the crossover temperature scales with the density of ex- 
cess quasiparticles. Another consequence of this nonequi- 
librium analysis is to understand the tunneling matrix 
elements of a superconducting qubit crucially depends 
on the coupling between the excess quasiparticle and en- 
vironment. This dependence even for low-lying quasi- 
particles is not negligible and that only for qubits with 
negligible superconducting phase fluctuations can be ap- 
proximated to the well-known sin 2 ip/2 matrix. This may 
modify, alter, or obviate the qubit transitions such that 
a forbidden transition in one qubit can become allowed 
in another one with different phase fluctuations. 

In a superconductor adding an electron-like excitation 
above the Fermi momentum (i.e. k > hp) or a hole- like 
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excitation below it (k < kp) will bring the superconduct- 
ing region out of equilibrium. A Gedanken experiment 
may help to understand why the region goes out of equi- 
librium |26j . Consider two identical superconducting is- 
lands placed near each other, each in parity odd state. 
Without loss of generality we consider each state carries 
a single excess quasiparticle. Once a contact between 
them is made, the two excess quasiparticles are recom- 
bined after exchanging phonons and relax to the conden- 
sate state such that in the overall island there remained 
no unpaired quasiparticle. This signals that the initial 
odd state is not a true thermodynamic equilibrium state 
because if so, after the contact is made the final system 
would have preferred to maintain the initial state of its 
subsystems. 

Let us now consider an electron (or hole) of mass m and 
momentum k has the kinetic energy ^ k = h 2 \k\ 2 /2m — (i s 
from the condensate chemical potential fi s . In supercon- 
ductors the quasiparticle excitations are positive energy 
levels at E k — + A 2 . Based on the BCS theory the 
occupation probability of the cooper pair state (— k, k) is 
v\ = (1 — £k/Ek)/2 and similarly its emptiness probabil- 
ity = 1 — v\. The total charge in a unit volume is 
q — 2e J2k u kfk + u fc(l — /-fe)i where e is the elementary 
charge unit, the sum is over momentum states k above 
and below kp, and the factor 2 indicates the spin degen- 
eracy. The function f k is the occupation probability of 
state k. This function in the thermal equilibrium is the 
Fermi-Dirac function /P. In the definition of q, the sum- 
mand represent the occupation probability of the elec- 
tron state k. The first term is the probability that in the 
lack of cooper pairs in the state (— k, k) a quasi-electron 
occupies k. The second term is the probability that in 
the presence of the pair state (—k,k) the quasiparticle 
state — k is empty of quasiparticle. The charge q can be 
rewritten in the form of a summation q = q s + q qp for 
the definitions of supcrfluid charge q s — 2e^ fc ?;^, and 
quasiparticle charge: 

<?qp = ^e^ulfk - v 2 k f- k . (1) 
k 

In equilibrium the energy distribution function (i.e. 
/j?) is an even function since the occupation probability 

for the states k and — k are equivalent; i.e. f k =f^ k - For 
such a function the effective quasiparticle charge can be 
simplified into q qp = 2J2 k Qkfk, where q k = e(u\ - v%) = 
e£ k /E k is the charge imbalance per unit volume. This 
imbalance is maximally +e at £ 3> A and by decreas- 
ing momentum it is suppressed to zero at kp, turns into 
negative values for k < kp, and finally reaches — e at 
£ k *C —A. A zero effective quasiparticle charge expected 
with even occupation probability functions, thus in equi- 
librium the total electronic charge is restricted to that 
of superfluid. This explaines why equilibrium quasiparti- 
cle tunneling are depleted from superconductors as they 



are electronically neutral. However, in nonequilibrium 
quasiparticles the energy distribution function (f k ) may 
carry an odd term whose contribution may give rise to 
a nonzero charge imbalance (q qp ^ 0). Let us indicate 
that the total charge is suppressed by Coulomb interac- 
tions, which implies that any change to the quasiparticle 
charge are accompanied by equal and opposite change 
in the pair charge |27) . An injected current in a tun- 
nel junction between a normal metal, held at bias volt- 
age V, and a nonequilibrium superconductor is given by 
the sum of the usual Giaever tunnel current I(V) and a 
bias-independent extra current that carries the contribu- 
tion of charge imbalance in the form of — qk{fk ~ f k )- 
This current is a physical quantity and therefore indi- 
cates that the charge imbalance nonequilibrium effects 
lead to a measurable steady-state voltage, (for details 
see [231 [24]). 

In quantum information processing devices the con- 
cerns about the noises forces to put the qubits at the 
temperature less than 1/10 of the superconducting crit- 
ical temperature. Parker in [28] showed that at such 
a low temperature the recombination time takes longer 
than the equilibrium thcrmalization. This recombina- 
tion time in superconducting qubits has been measured 
to be of the order of milliseconds, which is longer than the 
thermal relaxation time by at least two orders of mag- 
nitude. This poses a bottleneck causing a rise in the 
number of nonequilibrium quasiparticles in the supercon- 
ductor. Charge neutrality ensures that the charge imbal- 
ance induced by these quasiparticles is compensated by 
suppressing the superfluid charge q s . This is in fact the 
reason the charge imbalance does not dramatically affect 
the macroscopic quantum dynamics of the junction. The 
decrease of superfluid charge implies the suppression of 
chemical potential of superfluid from /i s to /i s — S^i s [29j . 
Consequently after the entry of a quasiparticle into su- 
perconducting leads and islands their chemical potentials 
decreases. 

For superconducting state under the influence of ex- 
ternal dynamic pair breaking, one can look at the un- 
paired quasiparticles as Fermi gas in equilibrium and ap- 
proximate its nonequilibrium occupation probability by 
/o (e, jU — Sfj) [Ml HHHH] • This function was previously 
used in the literature, for instance to determine the ex- 
perimental tunnelling rates of elastic and nonelastic ex- 
cess quasiparticles between two unentangled cooper-pair 
boxes [15 . Similarly, Palmer et.al. in [323 used this 
nonequilibrium function and described to a surprising 
accuracy in the differences between odd and even state 
plateaus in the Coulomb blockade staircase of a charge 
qubit. 

Let us consider a density of equilibrium and nonequi- 
librium in a region of superconductor. As discussed 
above, only the nonequilibrium density n nqp contributes 
to the chemical potential nonequilibrium. By definition, 
the number of nonequilibrium ones in a volume V can be 
written as N Dqp = V J deD(e)[f (e, /J, - 6/j) - /o(e,A*)], 
where D(E) is the density of states D(Ef)Q(\e\ — 
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A)|e|/-\A 2 — A 2 with step function &(x) — 1 for x > 
and otherwise. Substituting the Fermi-Dirac functions 
as well as the density of states in superconductors, one 
can find the chemical potential shift in the section i of a 
superconducting device is (see Appendix A for details): 



Sfii = fc^Tln 1 



1 + 3k B T/8Ai 



^i/k B T 



(2) 

One can see from eq. ^ that lowering temperature 
brings the chemical potential to further decrease. Once 
the chemical potential is shifted more than the qubit 
temperature, the quasiparticle tunneling enters nonequi- 
librium regime. From eq. ^ one can find this onset 
temperature is about T Q = 2A/fcs(lnz — ln(lnz)), where 
2 = 12n(AD(E f )/n nqp ) 2 . Note that n nqp /AD(E f ) is 
the ratio of excess quasiparticle per cooper pair in the 
superconducting region, which is fabrication-dependent. 
In a device with more excess quasiparticles tunneling 
the device enters the nonequilibrium regime at higher 
temperature. Let us immediately check the relevance 
of this prediction by plugging in typical parameters for 
qubits. In Ref. [TS] a quasiparticle density of 10//im 3 
was experimentally measured in charge qubits. For alu- 
minum superconductor that is typically used, we have 
A/k B = 2. IK and AD(E f ) = 2.8 x l0 6 /pm 3 . This yields 
the onset temperature at 0.16K, which is experimentally 
the temperature below which the tunneling rate becomes 
deviated from the exponential decay of equilibrium in ex- 
periment (see Fig. 2b in Ref. [17]). As a result, once the 
onset temperature is known it helps to estimate the den- 
sity of excess quasiparticles n nqp in a superconducting 
qubit. 

Nonequilibrium quasiparticles when tunnel in and out 
of a superconductor, their induced charge q qp may cou- 
ple them to the environment. As described in Fig. 
the change of parity in the superconductor from even to 
odd state yields decay or absorption of energy exchanged 
between qubit and the Fermi gas of quasiparticles. The 
environmental coupling of tunneling quasiparticles cou- 
pling induces fluctuations in the superconducting devices. 
In qubits with macroscopic superconducting phase ip this 
causes the phase to fluctuate (f>(t) = ip + S(p(t). The ef- 
fect of these fluctuations on the tunneling process has 
been studied in weak coupling [TTJ UH] and strong cou- 
pling limits |19| . The cooper pair phase fluctuation 8<j)(t) 
is defined through the noise SV(t) of voltage across the 
junction S<f>(t)/2 = (e/H) J* dt'SV(t'), where the phase in 
this paper is the conjugate to the Cooper pair charge q s . 

Using the standard methods in the environmental P(E) 
theory, the noise power of phase fluctuation, i.e. S(t) = 
(S(j>(t) 6(f)(0)) , can be determined from the fluctuation- 
dissipation theorem the way how it describes Johnson- 
Nyquist noise [33] . The quasiparticle channel is assumed 
to be in parallel to supercurrent, junction capacitance, 
and an external impedance channels. The noise from 
the resistor gets shunted through the Josephson junc- 



tion at low frequencies therefore the electromagnetic en- 
vironmental effects can be encoded in the transmission of 
quasiparticle channel. Considering a qubit junction with 
high quality factor, the effective impedance of the junc- 
tion is restricted to the absorption and emission of a sin- 
gle mode Up, i.e. Z e g = tt/2C[S(u> — uj p ) + S(uj+u!p)}. The 
junction plasma frequency is uj p — ^2itI c /CQq where 
I c is the critical current and C the papacitance of the 
junction and 4>o is the quantum of flux <E>o ~ 2.07 x 
10~ 15 Wb. Let us introduce the dimensionless parameter 
wave impedance p c = 4ttZq /Rk where Z — y/Q /2Tr I C C 
and the electronic resistance quantum Rk = h/e 2 — 
25.8kl7. Note that in this definition the cooper pair quan- 
tum resistance (i.e. Rk /4) was considered, which makes 
it 4 times bigger than p definition in Ref. [33J . We encode 
the difference in the subindex "c" . The vacuum fluctua- 
tions of the condensate is controlled by its corresponding 
wave impedance, i.e. S(0) = p c . The phase in our anal- 
ysis is conjugate to the cooper pair charge, thus in our 
derivation of tunneling rate from Fermi's golden rule |19j 
the phase fluctuations appears in the form of 5/4; (for 
details see Appendix B). From the above definitions one 
can drive p c = \J E c /2Ej, for the cooper pair charging 
energy E c , and the junction energy Ej. 

From Fermi's golden rule one can derive the tunneling 
rate for any strength of environmental effects and the 
broader energy range enforced by the multi-boson tran- 
sitions [19]. For a typical superconducting qubit with 
arbitrary coupling and bias phase, this results into the 
rate of eq. ( B8 ) in Appendix B. This can be simplified in 



the case of a phase qubit with large capacitance, where 
the phase is localized around its bias value at zero: 



r = r, 



7T 



(3) 

where we defined £ = Eif/2kgT and the cofactor r o := 
y/WfEjAEj n nqp /8D{E f )R N e 2 . The function K n (x) 
is the Modified Bessel function of the second kind. 

Neglecting nonequilibrium effect in the chemical po- 
tential makes the corresponding rate to scale with tem- 
perature as exp (— A/ksT) ^/ksT. At low temperature 
the equilibrium tunneling rate is heavily suppressed and 
effectively blocked. However, our nonequilibrium results 
indicate not only an enhanced rate at low temperature, 
but also a finite nonequilibrium quasiparticle relaxation 
time at zero temperature, which are consistent with ex- 
periments pH51 fT5HT7] . For a transmon, Ej/E c « 30 and 
for a flux qubit « 50 the correction from phase fluctua- 
tions are large enough to be visible in deformation of the 
quasiparticle rate function. 

Let us now consider a rather more general case of tun- 
neling nonequilibrium quasiparticles in a qubit of bias 
phase ip. In the limit of Eif kf,T, from from eq. (B8) 



one can determines the NQP tunneling rate in a typical 
superconductor qubit: 
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FIG. 1. (a) The transition between even and odd parity 
states; both initial and final states are tuned at a fixed 
charge gate n. The red arrows indicates the parity transi- 
tion in the qubit as a result of quasiparticle poisoning, (b) 
The energy profile at a junction between a large reservoir 
(lead) and a small island. The even to odd transition rate 
r eo and the opposite one are depicted by arrows from their 
corresponding initial to final states. Large transition energy 
makes the even to odd transitions much more frequent than 
opposite, (c) Temperature dependence of NQP relaxation 
time. The superconducting energy gap is A = 2.1 K, nor- 
mal junction resistance R = 10K57, the nonequilibrium quasi- 
particle density n nqp = 0.1/pm 3 , and the density of states 
D(Ef) = 10 6 //im 3 K. The junction is coupled to environment 
in the range of p — 0.02 to 0.2 with steps of 0.02 from top 
to bottom, respectively; here p := p c /4. The highest curve 
corresponds to Ej w 75E C and the lowest one to Ej w E c . 
Insets: (upper) Temperature dependence of NQP relaxation 
rate, from nonequilibrium analysis (solid) compared to that in 
equilibrium (dashed), where the upper curves is for p = 0.01 
and the lower ones for 0.4; (lower) the slope of NQP relax- 
ation rate as a function of temperature and the environment 
coupling p limits, Positive rate slope below the curve p* and 
negative slop above it. 



r = 



y/E c /32Ej + 1 / AE C 



^nqp 



8e 2 i? 



N 



E if Ej Di(E f ) 



(l—r/costp) 



(4) 

where <q = (1 - p c /4)/(l + p c /A) « 1 - y/E c /8Ej. This 
indicates that in typical superconducting qubits the tun- 
neling matrix clement of computational qubit is deter- 
mined by the operator 1 — 77 cos (p and is not separable 
from the quasiparticle contributions. Only in qubits with 
negligible phase fluctuations (i.e. Ej 3> E c ) one can 
retrieve the well-known sin 2 (( i 5/2) tunneling matrix ele- 
ments, similar to Refs. |llH13j ; however note that even 



in this case due to the consideration of nonequilibrium 
occupation probability in our analysis the quasiparticle 
amplitude is different. 

In a qubit of large capacitance and transition en- 
ergy (i.e. Eif k^T) the temperature dependence of 
nonequilibrium quasiparticle relaxation time can be di- 
rectly determined from eq. i.e. 




(5) 



At zero temperature this relaxation time shows to be 
finite and fabrication-dependent as it scales as T nqp ~ 
(l/n nqp y/Ei f E c /AEj)expy/E c /8Ej. From this one 
can expect in a device with poor shielding against the 
quasiparticle poisoning the relaxation time suppresses. 
Moreover, in a flux qubit of Ej/E c ~ 50 has a relaxation 
time almost three times longer than that in a quantro- 
nioum qubit of Ej/E c s» 5, and almost the double of 
that in a transmon, making qubits with larger E c /Ej 
less sensitive to the tunneling of excess quasiparticles. 
Note also that in qubits with smaller qubit transition en- 
ergy decreasing temperature makes the tunneling process 
remarkably slower. 

Figure ([lj;) displays temperature dependence of quasi- 
particle relaxation times from eq. ([3| in a qubit with 
large capacitance. The junction parameters indicated in 
the caption are for typical superconducting qubits. Here 
we assumed the allowed parity changing transition fre- 
quency is 1 GHz and the environmental phase fluctua- 
tions are in the range of Ej rj 75E c to Ej w E c in 
equal steps from top to bottom, respectively. By de- 
creasing temperature below the onset temperature at 
0.16K, the relaxation time gives up its exponential in- 
crease and turn into a regime with much slower increase. 
This precludes reaching the millisecond timescale ob- 
served in experiments. In qubits with weak coupling 
to environment p c the slope of relaxation time is neg- 
ative, making slower tunnelings by decreasing tempera- 
ture. However, in junctions with smaller capacitance and 
consequently larger phase fluctuations, the slope turns to 
positive after passing a maximum time. Both of these 
cases are consistent with experiment |16] 117] . There is 
a crossover environmental coupling p*(T) below which 
the relaxation time slope is positive and above which 
negative. At zero temperature the critical coupling is 
p*(Q ) — (A/6Eif — 1/8) ~ , (see Appendix C and eq. 
(CI) for finite temperature result). 

In conclusion, we developed a theory for the tunnel- 
ing of nonequilibrium quasiparticles across a supercon- 
ducting junction using the chemical potential disequili- 
bration that quasiparticle charge imbalance induces on 
superconductor. We calculated the transition rate for 
arbitrary bias phase and phase fluctuations and showed 
at low temperature the rate is enhances by some orders 
of magnitude compared to thermal equilibrium. The on- 
set temperature for nonequilbrium effects can be used 
to estimate the density of excess quasiparticles tunnel- 
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ing. Moreover, a finite energy relaxation time is pre- 
dicted from this model. Predictions of this model are 
in reasonable agreement based on a prior experimental 
measurement and lead to further analysis of the noncqui- 
librium quasiparticle tunnelings. Our model shows that 
the fluctuations of superconducting phase affect the qubit 
tunneling matrix elements to be 1 — rj cos ip and this ma- 
trix only in the limit of negligible fluctuations is simpli- 
fied to the well-known sin 2 (p/2. This can dramatically 
modify the allowance and forbiddances of transitions in 
the qubit during a tunneling. With a good understand- 
ing of the physics going on in a superconducting qubit, 



more experiments can be performed to further improve 
the coherence of superconducting qubits. 
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Appendix A: Chemical potential shift S/i 

Using the definition of the nonequilibrium quasiparti- 
cle density n nqp = N nqp /V = J deD(e) [f (e, fi - 5 ft) - 
/o(e,/x)] one can use the definition of density of states 
in superconductor based on the normal metal density of 
states D n and in the limit of low temperature finds 



E 



-Mr. 



l-X 



7 2A 



V-B 2 - A 2 

(! ~ X ) ,,-A/fcT 



X 

V2nAkT ^1 + 



l + X + e- E / kT + Xe E / kT 
( 1 J3_ 
A WE- A + 4A 
3kT\ (l-X) A/kT 
8A X 



dE 



D n {E f ) 



4e- 



i 

e (E-A)/kT 



dE 



D n (E f ) 



where we used the definition X — exp(5/i/kT). Solving 
the final equation results into the chemical potential shift 
ofeq. ^. 



Appendix B: Rate 

We formulate the mathematical description of quasi- 
particle tunneling in an environment. We start from the 
Hamiltonian 



H = H 



BCS,1 



H 



BCS.2 



H T + H C 



(Bl) 



for BCS Hamiltonians in mean field form describing the 
two electrodes and the tunneling Hamiltonian 

Ht = J2 ( T fc*4<r,ic^2 + T fe * c^a^i) (B2) 
hi 

By diagonalize the BCS Hamiltonians through the Bo- 
goliubov transformations, one can rewrite the tunneling 
Hamiltonian, eq. (B2) in the following form 



Bt= Yl ( T ki\wi 2 \e^ /2 -Tl^ k \v kl v l2 \e-^ 2 ) 

fe,CT=t4 

X (lt al %a2 + 7l feCT2 7-i-i) + H T 2 (B3) 



where we introduced the phase difference cj> — (pi — (f> 2 . 
The term Ht2 contains operators that change the num- 
ber of quasiparticles by two, i.e., contain terms of the 
structure 77 and -yt^t hence changing the total number 



of quasiparticles in the setup, which do not contribute to 
the quasiparticle rate -these terms contribute to Joseph- 
son and Andreev processes. 

We now want to evaluate the Fermi's golden rule rate 
for a transition that transfers a quasiparticle from the 
electrode 1 to the electrode 2. The relevant matrix ele- 
ment is 

((N, - l) ka , (N 2 + \H T \ N lk(T ,N 2la ) = 
T kl e 1 ^ 2 \u kl u l2 \ - Tl^e-' 1 ^ 2 \v kl v l2 \ 

In order to capture the influence of the environment, 
we apply the ideas of P(.E) -Theory [33]. There, the envi- 
ronmental Hamiltonian is described by an oscillator bath 
-ffbath = Y^n a -'n^n a "n an d couple the oscillators linearly 
to our quasiparticle system H env _ c = NieSV, where we 
assume that only the first reservoir fluctuates - this cor- 
responds to a specific choice of gauge. The total number 
and voltage operators are 



E 



-fco-l 



Cfco-1, 



sv 



A, 



(B4) 



In order to capture the influence of the environment, 
we apply the ideas of P(P)-Theory [34]. There, the envi- 
ronmental Hamiltonian is described by an oscillator bath 
-ffbath = Yin aj n^n a "n an d couple the oscillators linearly 
to our quasiparticle system H env - C — NieSV, where we 
assume that only the first reservoir fluctuates - this cor- 
responds to a specific choice of gauge. 

In the limit of single photon exchange in a large capac- 
itance junction, one can rewrite the quasiparticle contri- 
bution as 



e 2 R N 



e- 2 Pp 



dEh(E){l- f 2 {E + hu)) 



E(E + hio)- Ai A 2 cos ip + [E (E + hw) + Ax A 2 cos cp] S(t) /4 



\ 



{E 2 ~ A 2 ) [{E + tuo) 



A 2 

^2 



(B5) 



In quantum information processors the cutoff of the 
Fermi functions are well below the edge of the gap, so 
f (E- 6^)11- f(E+E if -5^)} « expI-^-^O/fcsT]. 
In Josephson junction with two reservoirs we assume the 
two chemical potential shifts are the same. If the junction 
is between a large lead and a small island the tunnel rates 
are dominated by the quasiparticle density in the leads, 



rather than in the island and here we assume the label 1 
is assigned to the larger lead. 

Taking the integrals of all terms and simplifying the 
rate after some steps of using calculus one can find the 
tunneling rate of nonequilibrium quasiparticle for arbi- 
trary superconducting phase difference ip across the junc- 
tion and parity transition energy Eif in a junction made 
of the same superconducting material: 



T eo [p,T,E lf ,S(t)} = -l^e^- 2p pAet-^ 



e 2 R N 



1 - cos<p kT 



3 + cos <p 9 — cos tp kT 



1 + cos ip 



K (0 + £ 



4 32 
kT / 3 - cos tp 9 + cos <p kT 
X v i 32 X 



A 



Ki(0 

S(t) 
4 



(B6) 



One limitation is to consider that the transition energy 
between the two parity states is of the order of the qubit 



energy. Therefore, in the limit of £ = Eif /2kT 1 given 
that Eif/A -c 1, from eq. (B8) one can derive: 



T eo [p,T,E lfl S(t)] 



^- 2 P, 



— e kT 
In 

3 — cos ip 



( ( 3 + cos ip 



e kT 



2 

kT^ 



4 

COS ip 



kT^l- 



1 — COS ip 



3 + cos ip\ kT 
A" 



3 /3 - cosp\ kT\ A 



A 

7! 



(B7) 



In the low temperature limit and large superconduct- magnetic vacuum fluctuations. In this case we can sub- 
ing gap, using eq ([2| one can consider the fluctuation of stitute further simplify the tunnelling rate of eq. (B7) 
phase is minimal and therefore it is of the orer of electro- into: 



r eo [ip,E if> T = 0] = 



AE C 
EifE., 



8e 2 R N Di(Ef) 




1 - 



E c 
2Ej 



cos p 



(B8) 



Appendix C: Crossover coupling p*(T) 

In order to find the condition for existence of the min- 
imum tunneling rate at low temperature let us take a 



partial derivative from the tunnelling rate of eq. ([3]) 
with respect to the temperature. Simplifying the equa- 
tion provides the corossover environmental fluctuations 
is 



7 



p*(T) = 



2A 

Eif 



£J*i(fl- 



I (#o(0 



A 



2A 

Si/ 



^2(0) 



(Cl) 



Note that the crossover p depends nontrivially on the 
superconducting phase. If one uses eq. ( B8 1 the phase 
dependence of the minimum rate can be found. 
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